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Abstract 

Several time-independent methods have been devised for measuring the phase 
7 of the Cabibbo-Kobayashi-Maskawa |||] unitarity triangle. It is shown that 



fT^ ■ such measurements generally suffer from discrete ambiguity which is at least 

> 

rn . 8-fold, not 4-fold as commonly stated. This has serious experimental implica- 

P^ ■ tions, which are explored in methods involving B -^ DK decays. The mea- 

o' 

0> ■ surement sensitivity and new physics discovery potential are estimated using 
a full Monte Carlo detector simulation with realistic background estimates. 
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I. DISCRETE AMBIGUITIES 

Direct CP-violation within the standard model takes place when two or more amplitudes 
with different CKM phases interfere in a time-independent manner. As an example, take 
the decay of a B'^ meson to a final state /"'", which can proceed through two amplitudes. 
The partial decay widths can be written as 

r(5± ^ /±) = Al + Al + 2A1A2 cos{S ± (/)), (1) 

where Ai and A2 are the magnitudes of the interfering amplitudes, is the CKM phase 
difference between the amplitudes, and 6 is the CP-conserving phase difference. CP-violation 
is often considered in terms of a non- vanishing CP-asymmetry, 

A _ ^(B+ ^ /+) - TjB- ^ /-) _ 2i?sin,5sin0 

T{B+ ^ f+)+T{B- ^ f-) l + R'^ + 2Rcos6cosf ^^ 

where R = A^j A^. It is clear that A gives a determination of 7 which has a 4- fold ambiguity, 
due to its invariance under the two symmetry operations 

"J exchange • 7 ^'^ "B- {^) 

These ambiguities were also noted in methods which do not rely on a decay rate asymme- 
tryi. 

The S'sign ambiguity is generally considered unphysical, since application of ^sign to 7 



values within the currently allowed range |T6 



g = {40° < 7 ^ 100°}, (4) 

yields values which are not within g, and are therefore inconsistent with the standard model. 
We note, however, a third symmetry of Equation (|l], which makes the ambiguity 8-fold: 

'S'tt : 7 ^ 7 + TT, 6b -^ Sb -n. (5) 



This symmetry has so far been overlooked, perhaps because it does not appear to be a 
symmetry of A 

The significance of Sj^ arises from the fact that while 3-,^ and ^sign may be unphysical 
within the standard model, the operation S-^Ssign, applied to 7 G (7, results in values which 
are in or close to g. Thus, given the finite sensitivity of future experiments, it may be 
impossible to deem S'7r5'sign unphysical, hampering the ability to test the standard model 
using such measurements of 7. 

In some 7 measurement methods, there is no a-priori knowledge of the magnitudes of 
some of the interfering amplitudes. These magnitudes are then obtained simultaneously wth 
7. In this case, additional "accidental" ambiguities will exist if more than one magnitude 
value is consistent with the data. An example is given below. 

II. QUANTITATIVE EXPERIMENTAL IMPLICATIONS 

B -^ DK decays, recently observed by CLEO 0, provide several ways to measure 
the weak phase 7 = Big{—VudV*^/VcdVX)- Since non-standard model effects are expected 
to be small in such decays, comparing these measurements with experiments which are 
more sensitive to new physics may be used to test the standard model [Q. Gronau and 
Wyler (GW) 0] have proposed to measure 7 in the interference between the h —>■ cus decay 
B^ -^ D^K^ and the color-suppressed, h -^ ucs decay B^ -^ D^K^ . Interference occurs 
when the D is observed as one of the CP-eigenstates D°2 — ~7^\^'^-^^i^ which are 
identified by their decay products. The interference amplitude is 



V2 A{B+ -^ DI2K+) = ^B{B+ -^ DOir+)e^(^^+^) ± ^JB{B+ -^ D^K+), (6) 

where 63 is a CP-conserving phase. The value of 7 is extracted from this triangle relation 
and its CP-conjugate, disregarding direct CP- violation in D^ decays 0. Several variations 
of the method have been developed f^J^l- 

In practice, measuring the branching fraction B{B^ -^ D^K^) requires that the D^ be 
identified in a hadronic final state, / = i^~7r+(r27r)°, since full reconstruction is impossible 



in semileptonic decays, resulting in unacceptably high background. Atwood, Dunietz and 
Soni (ADS) ||^ pointed out that the decay chain B^ -^ D^K^, D^ —>■ f results in the same 
final state as B^ -^ D^K^, D^ —>■ /, where the D^ undergoes doubly Cabibbo suppressed 
decay. Estimating the ratio between the interfering decay chains, one obtains 



A{B+ -^ D^K+) A{DO -^ f) 



Vcb Vus Qi 



^(^° - /) ^, 6 (7) 



The numerical value in Equation (^ was obtained using |V^b/l^Jfe| = 1/0.08 Q], |\4s/Ks| = 
0.22, lai/asi = 1/0.26 H, and 

^JD' - /) _ 0031 (S) 

which is the ratio measured for / = K^tt^ [l^. Equation (|^) implies that sizable interference 
makes it practically impossible to measure B{B^ -^ D^K^), and the GW method fails. 

ADS proposed to use the interference of Equation (0) to obtain 7 from the decay 
rate asymmetries in B^ -^ fiK~^, where /j, i = 1,2, are two D final states of the type 
i^^7r+(n7r)°. Measuring the four branching fractions, B{B^ —>■ fiK'^), B{B^ —>■ fiK~), one 
calculates the four unknowns B{B^ -^ D^K~^), 7, and the two CP-conserving phases asso- 
ciated with the two decay modes. B{B^ -^ D^K^) and the D^ decay branching fractions 
will have already been measured to high precision by the time the rare decays B^ -^ fK^ 
are observed. In addition to the similar magnitudes of the interfering amplitudes, large CP- 
conserving phases are known to occur in D decays [jl3|, making large decay rate asymmetries 
possible in this method. 

Jang and Ko (JK) [|l^] and Gronau and Rosner [|^ have developed a 7 measurement 
method similar to the GW method, but in which B{B^ -^ D^K"^) is not measured directly. 
Rather, it is essentially inferred by using the larger branching fractions of the decays B^ — > 
D^K^, B^ -^ D^K^ and B^ -^ Di2K^, solving in principle the problem presented by 
Equation (|^). 



III. COMBINING THE ADS AND THE GW METHODS 

Since the b -^ ucs amplitude m B ^ DK is very small and hard to detect, several 
methods will have to be combined in order to make best use of the limited data. Quantitative 
estimates of the resulting gain in sensitivity are rarely conducted, since they require realistic 
efficiency and background estimates, and depend on specific phase values. Here we undertake 
this task for the case of combining the ADS and GW methods (contributions of the JK 
method are commented on later). In this scheme, one obtains the unknown parameters 

^^[B{B+^D'K+), 7, 6b, 6n}, (9) 

where ^^i = aj:g[A{D^ -^ f)A{D^ -^ !')*]■, by minimizing the function 

,.2/^x _ /^ ^(0-Q>» V I ( aiO-am y ( biO-b„, y ( HQ-bm y , . 

^ ^^^ [ ^a.n, ) ~^[ Aa„ ) ^[ A6,„ ) ^ [ Ab^ J ^ ""^ 

with respect to the parameters ^. In Equation (0) we use the symbols 

am = B{B+ ^ fK+) 

bm = B{B+ ^ Dl,K+) (11) 

to denote the experimentally measured decay rates of interest, and 



B{B+ -^ D^K+) B{D^ ^ /) + ^B{B+ -^ D^K+) B{D^ -^ f)e 

1 , , 2 



JiSD+Ss+j) 



±JB{B+ -^ DOK+) + JB{B+ ^ DOir+)e*(^^+^) 



(12) 
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to denote the corresponding theoretical quantities, a^, b^, a(0 and 6(^) are the CP- 
conjugates of a^, bm, a(0 and 6(^), respectively. Axm represents the experimental error in 
the measurement of the quantity Xm- 

Several gains over the individual methods are immediately apparent: In the ADS method, 
a D decay mode is "wasted" on measuring the uninteresting CP-conserving phases. By 
contrast, when combining the methods, knowledge of a^, b^, am and bm in a single mode 
is in principle enough to determine the four unknowns, ^, even if Sd = Sb = 0. In practice. 



adding the D^2 modes will decrease the statistical error of the measurement. In both the 
GW and the ADS methods, the ability to resolve the S'exchange ambiguity depends on the 
degree to which 6b varies from one B^ decay mode to the other. Experimental limits on 



CP-conserving phases in i? ^ Dn, D*-k, Dp and D*p |T^ suggest that 5b niay be small, 
making the S'exchange resolution difficult. When combining the methods, however, we note 
that 6(^) and b{C,) are invariant under 7 ^^ 6b, whereas a(^) and a(^) are invariant under 
•y ^^ 6b + Sd- The ^exchange ambiguity is thus resolved in a single B~^ and D decay mode in 
which 6d is far enough from or tt. 

IV. SIGNAL AND BACKGROUND ESTIMATES 

We proceed to estimate the sensitivity of the 7 measurement combining the ADS and 
GW methods, at a future, symmetric e"*"e^ i?-factory, operating at the T(4S) resonance. 
The detector configuration is taken to be similar to that of CLEO-III ||T8[. The integrated 
luminosity is 600 fb~^, corresponding to three years of running at the full luminosity of 
3 X lO^'' cm~^ s^^ M with an effective duty factor of 20%. 



Crucial to evaluating the measurement sensitivity is a reasonably realistic estimate of 
the background rate in the measurement of a^, whose statistical error dominates the 7 
measurement error, A7. We estimated the background by applying reconstruction criteria 
to Monte Carlo events generated using the full, GEANT-based ||2^ CLEO-II detector sim- 
ulation. The event sample consisted of about 19 x 10^ e+e^ -^ BB events and 14 x 10^ 
continuum e^e^ -^ qq events, where q stands for a non-6 quark. Since the full simulation 
did not include a silicon vertex detector or Cerenkov particle identification system, these 
systems were simulated using simple Gaussian smearing. The Cerenkov detector was taken 
to cover the polar region | cos^| < 0.71. 

D^ candidates (reference to the charge conjugate modes is implied) were reconstructed in 
the final states K'lr^, K^ti^ti'^, and K^ti^ti^ti^. The 7r° and -D" candidate invariant masses 
were required to be within 2.5 standard deviations (cr) of their nominal values. A Dalitz 
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plot cut was applied in the K tt+tt" mode to suppress combinatoric background. The B^ 
candidate energy was required to be within 2.5o" of the beam energy. The beam-constrained 



mass, \J El — Pg, where Ei, is the beam energy and Pb is the momentum of the B^ candidate, 
was required to be within 2.5cr of the nominal B^ mass. Since the K^ and the D^ fly back- 
to-back, all charged daughters of the B^ candidate were required to be consistent with 
originating from the same vertex point. Continuum background was suppressed by applying 
cuts on the cosine of the angle between the the sphericity axis of the B^ candidate and that 
of the rest of the event, and on the output of a Fischer discriminant [H. In background 
events, the reconstructed K^ and K~ come from two different D mesons, or are due to ss 
popping, while signal events often contain a third kaon, originating from the other B meson 
in the event. As a result, 90% of the background events are rejected by requiring that an 
additional K~ or Ks be found in the event and be inconsistent with originating from the 
B^ candidate vertex. 

With the above event selection criteria, we find that continuum events account for over 
80% of the remaining background, with a rate of 7 events per 10^ charged B mesons pro- 
duced. This is comparable to the expected signal yield. Under such low signal, high back- 
ground conditions, significant improvement is obtained by conducting a multi-variable max- 
imum likelihood fit. In this technique, cuts on the continuous variables are greatly loosened, 
and the separation of signal from background is achieved by use of a probability density func- 
tion, which describes the distribution of the data in these variables. As has been the case 
in several CLEO analyses of rare B decays, we assume that the effective background level 



in the likelihood analysis, 5, as inferred from the signal statistical error, AS = yS + B, 
will be similar to the level obtained with the Monte Carlo simulation. Signal efficiency will 
increase, however, due to the looser selection criteria. 
The expected number of B^ —>■ fK^ signal events is 

Na = NB+a{i)e{K+f), (13) 

where Nb+ is the number of B^ mesons produced, and e{K'^f) is the probability that the 



final state be detected and pass the loosened selection criteria of the likelihood analysis. For 
given values oiSo, ^b and 7, we calculate a(^) using the D'^ — > K~t{^ , K~7t^it^, K^-k^-k^-k^ 
branching fractions from |12[, Equation (D, B{B+ -^ D^K+) = 2.57 x 10"^ 0, and B{B+ -^ 
D^K+) = 2.3 X 10"^ (obtained from B{B+ -^ D^K+) and the values used in Equation (0)). 

To estimate the efficiency e{K^f), we start with the values in |0, 44% for the K'tt'^ 
mode, 17% for the K^ti^-k^ mode, and 22% for the K^-k^ti^'k^ mode. These are multiplied 
by the efficiency of finding the third kaon (45%), and the particle-ID efficiency (68%). 
The particle-ID efficiency is composed of the probability that a well-reconstructed K^ be 
in the particle-ID system's fiducial region (83%), and that half the K~ daughters of the 
D meson also be in the fiducial region. The momentum of the other half allows good 
identification using specific ionization, as does the momentum of the third kaon in most 
events. An additional efficiency loss of 10% is assumed due to non-Gaussian tails, Cerenkov 
ring overlaps, etc. The final efficiencies are 13% for the K~7r+ mode, 5% for the K^-k^-k^ 
mode, and 7% for the K^-k^-k^-k^ mode. 

Since hm '^ dm, suppression and accurate knowledge of the background in the measure- 
ment of bfn is much less critical. Starting from the continuum background level in and 
applying vertex and particle-ID criteria, we arrive at a rate of 60 background events per 10^ 
charged B mesons. The number of signal events observed in this channel is 

Nt = Nb+ 6(0 eiK+) Y^BiD^ ^ q) e{c,), (14) 

i 

where e{K~^) is the efficiency for detecting the K~^ with the particle-ID criteria described 
above, and q are CP-eigenstate decay products of Di2- Using Table |, we obtain J2i B{D^ — >■ 
d) e(c,) = 0.011. 

V. MEASUREMENT SENSITIVITY 

To estimate the measurement sensitivity for given values of the "true" parameters ^ = ^°, 
we compute the average numbers of observed signal events using Equations (0) and (|I1|) . 
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An integrated luminosity of 600 fb ^ yields Nb+ = 640 x 10^. We assume that statistics 
will effectively triple if, in addition to B^ -^ D^K~^, one uses the modes B^ -^ D^K*'^, 
B+ -^ D*^K+, B+ -^ D*°K*+, 5° -^ D°K*° and fiO -^ D*°K*^. We therefore take 
Nb+ = 1900 X 10®. The resulting A'^^, A''^ and their CP-conjugates determine the experimen- 
tal quantities a^, bm, ^m and bm in the average experiment, ie., the experiment in which 
statistical fluctuations vanish. The minimization package MINUIT [^ is then used to find 
the parameters ^, for which x^iO is minimal in this experiment. Since the measurement is 
expected to be statistics-limited, only statistical errors are used to evaluate x^(0- 

To demonstrate ambiguities, the trial value of 7 is stepped between —180° and 180°, and 
Sd, 6b and B{B^ -^ D^K^) are varied by MINUIT so as to minimize x^(0- Such 7 scans 
are shown in Figure |l| for cases of particular interest. Evident from these scans is the fact 
that a large d'^x^iO/'^l'^ at the input value 7 = 7*^ does not guarantee that x^iO ^iii obtain 
large values before dipping into a nearby ambiguity point. As a result, the quantity that 
meaningfully represents the measurement sensitivity is not A7, but fcxc, the fraction of g 
which is excluded by the B —^ DK measurement, ie., for which x^(0 > 10- The larger the 
value of /exc, the greater the a-priori likelihood that predictions of 7 based on new physics- 
sensitive experiments will be inconsistent with the B -^ DK measurement, leading to the 
detection of new physics. 

To evaluate /exc, 540 Monte Carlo experiments were generated, using randomly selected 
input values in the range 7° G g, -180° < 5% < 180°, -180° < 5% < 180° (Note that in 
reality, the CP-conserving phases will be different in the different decay modes). Depending 
on the input phases, the numbers of observed signal events varied between 700 < Nb < 1050, 
< A'a < 130. For each set of phases, a 7 scan was conducted in the range 7 G (7, and /exc 
was taken to be the fraction of the area of the scan for which x^iO > 10- 

The /exc distribution of the 540 random experiments is shown in Figure |[ Also shown 
is the distribution of the 91 experiments for which | sin((5B)| < 0.25. /exc tends to be larger 
in this case, since small values of x^(0 associated with the S'exchange ambiguity (even if the 



ambiguity is resolved) are pushed away from the center of g. Since the distributions of 
phases used in the Monte Carlo experiments cannot be expected to represent the actual 
phases in nature, it is not meaningful to study the /exc distribution in detail. Nevertheless, 
Figure indicates that this measurement may reduce the allowed region of 7 by as much as 
70%. 

VI. DISCUSSION AND CONCLUSIONS 

We have studied in detail the measurement of 7 using B -^ DK at a symmetric B fac- 
tory. Use of this measurement to detect new physics effects is complicated by low statistics 
and an ambiguity which is at least 8-fold, not 4-fold as often stated. We show that com- 
bining the ADS and GW methods helps resolve the S'exchange ambiguity and decreases the 
statistical error, compared with the ADS method alone. The ambiguities associated with 
the S'sign and 5,^ symmetries are irremovable in measurements of this kind. Even when the 
^exchange ambiguity is in principle resolved, in practice it still deteriorates the measurement 
by reducing x^iO ("^^ other experimental quantity of significance). 

Being ambiguity-dominated, the sensitivity of future experiments should be evaluated in 
terms of the exclusion fraction /exc, rather than the weak phase error A7. With a luminosity 
of 600 fb~^, we find that the B — > DK measurement can exclude up to about /exc ~ 0.6 of 
the currently- allowed range of 7. 

With 3 X 10*^ B mesons, 100% efficiency and no background, JK find A7 in their method 
to be between about 5° and 30° for 40° < 7 < 100°. Using more realistic estimates and 
noting out comments above, one would conclude that combining their method with the 
ADS and GW methods, while probably useful for the actual experiment, will not result in 
a dramatic change in the predictions of our analysis. 
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0.0036 0.087 0.0003 

0.0061 0.28 0.0017 
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total 0.011 
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a) 61=4.9°, 6b=16.9°, 7°= 122 



b) 60=32.3°, 8l= 1 6.9°, j°= 1 22° 



-c) 60=32.3°, 6b=16.9°, 7°=86.6° 



d) 6d=156°, 63=40.1°, 7°=61.3° 
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FIG. 1. x^iO ^s ^ function of 7 for different values of the actual phases, 6^, 5%, j^. For each 
value of 7, x^(0 is minimized with respect to B{B^ — > D^K~^), 6d and 6b- The points 7 = 7^ 
and 7 = (^^ are shown by a solid and a dotted line, respectively. Some asymmetry and noise are 
due to the dependence of the fit on the initial ^ values, a) The 8-fold ambiguity of Equations @ 
and (|5|) is demonstrated for small 6^. b) Increasing 6d, the ^exchange ambiguity is resolved, c) 
With 7 close to 90°, the 5^ and Ssign ambiguities overlap, d) The ^exchange ambiguity is resolved, 
but an accidental ambiguity shows up at 7 ~ 28°, with B[B^ — > D^K^) at approximately 4/3 its 
input value. 
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FIG. 2. The /exc distribution of all Monte Carlo experiments conducted, and experiments with 
sin((5B)| < 0.25. 
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